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PAFEC75 ON THE A.U.W.E. ICL 19045* 
\ A USERS GUID 


\/ PRECIS 
¥ — 


1. PAFEC75 is the latest development in the PAFHC suite which carries out 
stress or thermal analysis on any structure by means of finite element methods. 
A ,uide to the method of use on the 190480 at A.U.W.E. is given. It is intended 
that later a dedicated minicomputer will become available. 


INTRODUCTION 


2. PAFHO (Program for Automatic Finite Element Calculations) is a suite of 
programs which is capable of oe out a large range of engineering calcula= 
tions on engineering structures Its development started at Nottingham University 
in the middle 1960's and it was Ona issued commercially in 1970 as PAFHC70. In 
late 1972 this was superseded by PAFHC70* and extensions were issued, when avail-— 
able. It was implemented at A.U.W.E. in 1976. (Ref. 1). 


3. In 1976 the PAFEC group at Nottingham formed a company (PAFEC Ltd) to exploit 
the system. They then produced PAFHC75 based on a revised layout for input data, 
designed largely by users in 1974, together with internal changes to the scheme. 
This version was implemented at A.U.W.E. during 1977 and was finally made avail- 
able to general users of the 1904S* in January 1978. It is intended that PAFHC75 
will eventually be mounted on a dedicated PAF EC /GRAPHICS computer, hopefully at 
the end of 1978. 


4, Whilst PAFHC75 is well documented (Refs. 2, 3 and 4) it is the variation 
implemented on the 19045* that is described as far as the user is concerned. The 
method if implementation is described elsewhere (Ref. 5). It is also assumed 
that the user is already able to use the computer. (Ref. 6). 

THE SYSTEM 
Be PAFHC75 uses a series of 10 phases each loading a program. The data is in 


modular form and if some parts are omitted default values are used. The phases 
are described in the table below. 


ranee anes Detailed description 
ls . description 
Data modules are read in, default values are inserted 
and the modules are placed onto backing store. The 
NODES module is expanded so that all mid=side nodes are 
included. 
PAFBLOCKs Any PAFBLOCK data is replaced by the full nodal 
coordinate and topological description of the complete 
IN. DRAW 
structure 


The structure itself is drawn. At this stage it is not 
possible to show any results such as displacements, 

stresses or temperatures since these have not yet been 
evaluated. 


mesh of elements. 
Pre=solution 


housekeeping 


In this PHASE the constraints on the problem are con= 
sidered and a numbering system for the degrees of free- 
dom is derived. 


IN.DRAW 
constraints 


This PHASE is very similar to the PHASE 3 except of 
course the constraints which have been applied are 
shown. Conversely the degrees of freedom can be indica= 
ted on a drawing. 


The stiffness (or other such as conductivity, mass 
matrices of all the elements are found and put onto 
backing store. 


Elements 


4. 


Table (Cont'd) 


- ae Detailed description 
description 
S The system equations are solved for displacements, 
temperatures or whatever happens to be the primary 
unknowns in the problem being tackled. 


OUT. DRAW The primary unknowns in the problem (i.e. displacements 
displacements | or temperatures) are drawn. 


The stresses are found. 


stresses 


8 
9 
OUT.DRAW Stress contour, stress vector plate etc. are produced. 


6. In the system at A.U.W.E., programs have been compiled in advance and use 
made of a facility for increasing the core requirement of a program whilst 
ruming in order to enlarge 2n array BASE (which is used by most routines and 
whose size depends on the problem). Due to the limitation of 150K core size, 
imposed by the installation manager, the greatest possible size of BASE (at its 
maximum in phase 7) is 55000. It is hoped that the dedicated computer will have 
a paged environment thus allowing larger sizes of BASE. This system of pre= 
compiled programs reduces the time taken for each job. It does however assume 
that certain element combinations are excluded. 


Te The elements are taken from a library (Refs. 35 4) and are mainly broken 
up into 5 categories: 


Be Beam elements. 

b. Plane stress, plane strain and axisymmetric eiements (2D). 
ce Plate bending and flat faced shell elements (SHELL). 

d. Three dimensional solid elements (3D). 

e. Temperature or Laplacian elements (THERMAL). 


The only combination assumed is between beam and shell elements which have a 
combined program. Any other combixation will require a program to be specially 
written together with a macro to run it. This can be done by JCA1 staff at 
A.U.W.E. if required. 


on PAFBLOCKs (Phase 2) is a system whereby a mesh of finite elements can be 
created by specifying the topology of a block and the layout of the mesh in the 
block (Refs. 3, 4). 


9. Using the series of phases enables the job to be stopped at certain stages. 
These usually coincide with a plotting phase so that the drawing can be inspected 
for any curious element distortions etc. It is also possible to run phases 6 

and 7 on their own as they are those which involve the most computer time and 
core size in their execution and this enables the longer jobs to run. 


10. At the time of writing this guide it is EASIDATA only (Ref. 3) which is 
implemented. This should be sufficient for most users but it is hoped that the 
full DATA system (Ref. 4) will be implemented later. Whether this will be on 
the 1904S* or the dedicated computer will depend on the time scale and demand. 


PREPARATION OF DATA 


41. The written programs preclude the USE. and RHAD.FROM. facilities 

(Refs. 3, 4). The USE. facility would require a specially written program. The 
READ.FROM. facility is catered for in one of the parameters to the macro to run 
the system. Also the control of the phasing is done by parameters to the 
running macro and not by means of PHASE and SAVE commands in the CONTROL module 
(Refs. 3, 4). 


12. Apart from these limitations the data is prepared in the same format as 
shown in the manuals (Refs. 3, 4). There are two points to make however. I+ is 
preferable that the CONTROL module is the first in the data and the FAST. READ 
facility does not give such a benefit in speed of reading and will not allow 
continuation lines. 


RUNNING PAFECT75 


13. Three macros are used in the system. These are RUN75 and CLEAR75 which are 
only used for PAFHC75 and PLOTMAC which is used in common with PAFHC70+ (Ref. 1). 


14. RUN75, as might be expected, runs the system and has the following para- 
meters: 


DATA <filename>: Where €filename> is the data file. Note that the 
absolute name is required if the file is not in <account> 
specified by DY.below. This parameter is only requgred - 
for Phase 1 and is ignored in other phases. 


Where <ident> is a string of up to 8 characters on which 
the output file names are to be based. If this parameter 
is absent the job identifier (or the first 8 characters 
if more than 8) is used. 


OUT «<ident> 


DY <account> °¢: Where account is that in which the files are to be 
created. The account must exist and be available to the 
user. This parameter is not needed when running in the 
main account. It facilitates the erasure of files later 
(but see CLEART5 below). 


sT™mn $ Where n is the phase at which the job is to start. If 
omitted the job tries to start at Phase 1. 

Pin : Where n is the number of the plot phase required. Such a 
parameter is required for each plot phase required unless 
STn above denotes a plot phase. 

ENDn : Where n is the last phase required. 

BLOCK : When PAFBLOCKS (Phase 2) are being used. Not needed if 


ST2 is used. 


BOH 
THERMAL 


}: 


One of these is required to indicate which type of clement is 
being used in Phases 6 and 9. If omitted 2D is assumed. 


DYNAMIC ) 

STATIC ) 

THERMAL): One of these is required to indicate the type of solution 
(Phase 7). Of course for thermal work the parameter is the same 
as for Phases 6 and 9 so only one is required! If omitted STATIC 
is assumed. They can 4il be contracted to the first three charac— 
ters. (i.e. DYN, STA, THE). 


Where n is the total job time required. By default 10 mimtes 
(the system default). 


Where n is the program time required for each phase. By default 
300 seconds (5 mins). This time is carried into PLOTMAC if used. 


BASHn Where n is the size of base required. By default 5000 is assumed. 


C LEAR This sets up a call to CLEAR75 (para. 17) which clears out the 
files created, except for those output for graph plotting. This 
can be shortened to CL. 


Saves the magnetic tape file(s) carrying information between 
Phases 6 and 7 which are normally erased. (Para. 24). SV can be 
Causes the large number of output files to be listed on the line 
printer. Normally only those for Phases 7 and 9 (results) or the 
current file at any failure are listed. LI can be used. 


Suppresses the issuing of ENDJOB in a successful background job. 
Ignored in the event of any error. (It is assumed that the rest 
of the job depends on a successful runt!). 


RETAIN Retains the monitor file in the event of an ENDJOB issued by the 
macroe (See para. 26). RET or RT can be used. 


There are only two restrictions in the order of the above parameters: 

a. The DY parameter should precede DYNAMIC (or DYN) (if present). 

bs The ST parameter should precede STATIC (or STA) (if present). 
Apart from these restrictions the parameters can be in any order. 
jo. The formidable list of parameters should not frighten users! They are not 
all necessary at oncee In fact for a static, 2D analysis run in the main account 
with BASE less than 5000, job time less than ten mimtes (and presumably progtime 
less five mimtes), no plotting, no pafblocks and no retained monitor file 


the command needs to bes 


RJ <ident>, :<account>, RUN75, PARAM(DATA<filename>) 


However if more than 5 or 6 parameters are needed it is preferable, when adding 
a job to the operators" background queue, to set up the call in a small file ie. 


IN FILENAME 
RUN75 DATA PAFDAT,OUT PAFANS, BLOCK,PL3, END3,BASE 10000,RT,JT 900,DY < account> 
RHE 
Then 
RJ <ideni>, <account>, FILENAME (where«<account> here is the main account) 


This facilitates the task of the operators when issuing jobs from their background 
queue. 


47. CLEAR75 is responsible for dealing with the debris left by RUN75. It is 
able to cope with the remains of more than one job and requires the following 
parameters: 


/eident> : The identifier used in creating the files. 


DY ce account>: Causes the erasure to take place in the appropriate psuedo 
account. 


PLOT : Causes the erasure of associated redundant plot files. These 
are not usually erased, in case of plotter/computer failure. 


15. The first two parameters can be repeated as necessary. If all the files 
are in one pseudo account only one DY is needed. Otherwise one is needed paired 
with each /parameter. For example. 


CLEART5 /UBFSBO1, DY: XYZA, DY: XYZB, /UBFSBO2 
ral) 


es based on UBFSBO1 in :XYZA 
and UBFSBO2 in :XYZB 


would clear all fi 


CLEAR75 /UBFSBO1, DY: XYZA, DY: XYZB, /UBFSBO2, /UBFSBO3 
would clear all files based on UBFSBO1 in :XYZA 
UBFSBO2 in :XYZB 
and UBFSBO3 in :XYZB 


19. FLOTMAC is called from inside RUN75 but can also be called by the user. Its 
use is described elsewhere (Ref. 1). There is one modification that has now beer 
made, no real mag tapes are used and a series of pseudo mag tape files is set up. 
This has no direct effect on the user but makes the life of the computer staff a 

lot easier!t 


20. One modification made at A.U.W.E. is to check on the maximum size of BASE 
used (as opposed to that asked fort). This is noted at the end of the output for 
each phase and listed to the monitor file (see para. 29 and Appendix A pp1 and 9) 
and also affords a check on the size of BASE obtained (which is not always the 
same as that asked for due to limitations on core size). 


21. In the event of any program failure the macro RUN75 is ended (a background 
job is terminated even if NEJ is present) and the current output file (if one 
exists) is listed on the line printer. 


FILES CREATED 


22, TReIdent> is the line printer output file. It has a generation number corres 


pondins to the phase and a language code ANSA. If the first file for any phase 
becomes full (either due to the limit of 8192 transfers or physical size) a second 
file is assigned with code ANSB — and so on up to ANSZ. 


23. Beicent>S is the disc file containing backing store information (referred to 
1s BS in literature). It is 240 kwords long. After each phase is run Beident>S 
copied into H¢ident>S. In the event of a successful run Heident>S is erased 
-u* in the event of failure (usually due to BASE being too small or the time being 
insufficient) it is left in the filestore to be copied back by the macro before 
‘he phase is rerun. This has been found to be necessary as a restart may find 
S<ident>S changedtt 
21, M<ident>T is a pseudo mag tape file to carry information between phases 6 
and 7. If this file is nearly filled, a second is assigned from inside the pro- 
gram and so on. The language code of the first file is FIRS, subsequent codes are 
AAMT~AZMT—BAMT=BZMT—CAMT and so on to ZZMT thus catering for 677 files. Usually 
erased after Phase 7 is run but see SAVE parameter for RUN75 (para. 14). 


25. “ident>Qn is the card image file output by each plot phase to be presented 
to PLOTMAC where n denotes the phase i.e. Z«ident>Q3 for Phase 3. It is the file 
from Phase 10 which sets the limit of 8 characters on <ident> (from the limit of 
12 characters for a filename). These files can have more than one generation 
runber due to the card image file being filled. 


25. Meident>N is the retained monitor file name. 


It should be remembered that <ident> is normally the job identifier unless 
OUT is being used and the files are all created in the pseudo account (if any) 
in which the job is run. 


EXAMPLE OF RUN 


23. The monitor file of the run and the listings of the output from the phases 
are to be found in Appendix A. The monitor file (pp A1—A3) shows the contents of 
:AAHC before the run (p A1), the effect of NEJ, the contents of :AAHC after the 
min (p A3) noting that the job is now in sAAHC. Then the effect of CLEAR75 can 
he seen when the LD command shows that the only files left are the data file and 
plotter output file. A further use of CLEAR75 with the PLOT parameter clears the 
plot file as well. RUN75 without any ST and DATA parameters gives rise to the 
DISPLAY and the job is terminated. Another RUN75 with DATAORB was included in 
the macro which would have given rise to the DISPLAY DATA FILE DOES NOT EXIST. 


29. It can also be seen from the monitor file that the maximum size of BASE used 
was 2990 in Phase 9 (p A3) followed by 2022 in Phase 8 (p A2), 1851 in Phase 6 

(p A2) and 1759 in Phase 7 (p A2). From page A1 it can be seen that the size 
required in Phase 1 is only 131. 


30. The plot of the displaced shape from Phase 8 can be seen in Fig. 1. From 
the output from Phase 8 (p A31) the SCALE OF DISPLACEMENTS should be 0.994E=5 
per om. There is also some sorting of the layout of the labelling in the bottom 
right hand corner which will have to be done. 


37. It should be noted that without the LI parameter the only listings (apart 
from the monitor file) would be from Phase 7 (pp A22-A27) and Phase 9 (pp A33 - 
AAJ) a total of 17 pages instead of 40. 
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